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I. INTRODUCTION therefore their complete analytical solution using known 
The analysis of asymmetrical transient processes of tabulated functions is practically impossible. 

synchronous machines becomes of great importance, An analytical study of asymmetrical transient 
especially for highly used powerful turbo and hydro processes of a synchronous machine operating on an 
generators. Аз a rule, modern heavy-duty generators infinite power network through reactance was carried out 
operate in fairly large power systems and therefore it is quite fully only in (Goldberg, Bul, Sviridenko, & 
necessary to analyze asymmetrical processes not of single- Helemskaya, 2001). The author provides analytical 
phase generators (to which a large number of works are expressions for the stator current only for a two-phase 
devoted and this issue is sufficiently covered in the press), short circuit of the generator and a single-phase short 
but the operation of generators that lose connection with circuit to the neutral of the system and the generator. 


the power system (for example, with single-phase or two- 
phase short circuits, with phase-by-phase automatic 
restarts, etc.) (Basmanov & Votinsev, 2021) (Biswas, Pal, 
Werho, & Vittal, 2021) (Dong, Tian, & Ding, 2021). 
Additional difficulties for the analysis are presented by 
taking into account the reactivity connected between the 


The difficulties of solving a complete system of 
differential equations with periodic coefficients are 
overcome by numerically solving the problem using 
computer technology (CT). This work is devoted to the 
creation of a technique for numerical analysis of some 
transient asymmetric modes of a synchronous generator 
operating on an infinite power network using CT (Hugo, 


and power transmission. dde (ETL)):. Differential Gloria, Alvaro, Jesus, & Nikolas, 2022) (Jaramillo Serna 
equations describing transient asymmetric processes of a 


synchronous machine contain periodic coefficients and 


generator and the system (the reactivity of transformers 
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& López-Lezama, 2019) (Krishna, Sasikala, & Ganesh, 
2017) (Pai, 2014). 


The article presents the basic equations for 
calculating transient asymmetric modes on a CT, 
information on their programming, as well as the results of 
calculations on a CT. 


Il. BASIC EQUATIONS FOR CALCULATING 
ASYMMETRICAL MODES OF A 
SYNCHRONOUS GENERATOR 


The initial equations for the stator winding are 
compiled in the system of axes a, b, c; for rotor quantities, 
the system of axes d, q, O is used. Periodic coefficients in 
the equations of a synchronous machine are calculated at 
each interval of the numerical solution of the equations as 
a function of the angle between the stator and rotor axes. 
In this case, its periodic coefficients are expressed only in 
terms of and calculating them on a CT is not difficult. 
(Dong, Tian, & Ding, 2021) (Grigsby, 2001) (Guo, Bao, 
Xiao, & Chen, 2021) (Guseinov & Ibrahimov, 2012) 


The advantage of the used axis system over any other 
is that all current values correspond to real values and do 
not require recalculation of the results to obtain phase 
values. 


When drawing up equations for calculating 
asymmetrical transient modes of a synchronous machine, 
the following assumptions were made: 


1. The phase windings of a synchronous machine are 
symmetrical, i.e. they have the same number of turns, 
active resistance and mutual shift of magnetic axes. 


2. When considering the magnetic fields of the self- 
induction of the stator windings and the mutual induction 
of these windings with the rotor windings, only one 
harmonic of this distribution is taken into account. 


3. The magnetic permeability of the machine's 
magnetic core steel is assumed to be constant. Saturation is 
taken into account by choosing constant saturated 
parameter values. 


4. It is accepted that on the rotor, in addition to the 
excitation winding, there is one damper circuit along the 
longitudinal and transverse axes. 


The initial equations for calculating the modes of a 
synchronous machine in coordinates a, b, c are the 
following differential equations for the stator winding 


voltages: 
PY, = €, =, 
py,-e,-ir t. (1) 
PY. = €. -Lr. 
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where w,.w,.v,-— flux linkage of stator winding 


phases;  i,,i,,i,— stator winding phase currents; 
r,,lj,l,— active phase resistance of the stator winding; 


€,,€,,€,— Voltage at the terminals of the generator stator 
ee d : МУ : 
phase windings; p — "E differentiation operator with 
T 
respect to synchronous time т = 27ft . 


To this system of equations one should add the stress 
equations for the rotor circuits and the rotor motion 
equations: 


ПИР; 
PV ra = hea ka 
pU hd. $, OD 
pS =+ (M, M.) 
H 
po=S 


V (V uas к 

where 1,1, [— current flux linkage and active 
Гу, Туа» ы 

resistance of the excitation winding and damper circuits 
along the longitudinal and transverse axes; е; – voltage 
applied to the excitation winding; $— slip; H — inertial 
constant in electric rads; M,,— torque of the load on the 
shaft of a synchronous machine; M,- electromagnetic 
torque of synchronous machine; 0— working angle (the 
angle between the transverse axis of the rotor and the 
representing vector of phase voltages). 

To solve systems of equations (1) and (2) using any of 
the well-known numerical methods of Runge-Kutta, 
Adams Euler, etc. (Zakaryukin & Kryukov, Complex 
asymmetric modes of electrical systems, 2005), it is 
necessary that the number of variables equals the number 
of equations. Experience shows that it is expedient to 
express all currents through the flux linkages of the 
circuits. For this purpose, known relations are used, 
obtained in the calculations of symmetrical modes using 
the Park-Gorev equations (Goldberg, Bul, Sviridenko, & 
Helemskaya, 2001) (Lupkin, 1985) (Kryuchkov & others, 
2009) (Kundur) (Yusifbeyli, 2019). 


i, - ay, - by, -cVya 
i, -gV,-hy,, 

i, by, tdy, - ev, 

i, = суа - ev, fv, p (3) 
ia 7 hw, tky,, 
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where coefficients a,b,c,d,e,f,g,h,k are expressed 


through machine parameters as follows : 


2 2 

pa Se Q X ur S = Xu D 

^ А XX, -Xa 
b Xa & ka 7 Xu _ X Xu Ka h= Xaq 

А” А 7 X Xy x 

Xa X у Хаа Х„Х, Хаа X, 
A” A? XXu ta. 
(4) 
2 
A-X Aba Xa — х2.) Xaa (ea Xa TAS 


2 
— Хаа (хы X, — Хаа ) 
The parameters included in these expressions represent 
the mutual or total reactivity of the circuits: 
X, 
Xa = Ха tX; X, = Xag Xp 


=X Xf) Ха =Xaa t Xka) X rg = Xaq + Ху; 


To transition from stator currents i,,i,,iọ to phase 
quantities i,,i,,i,, we use the known relationships 


equations (Goldberg, Bul, Sviridenko, & Helemskaya, 
2001) (Lupkin, 1985) (Prabha & Lei) (Soldatov & Popov, 
2004) (Soldatov & Popov, 2005) (Zakaryukin, Kryukov, & 
Le, 2013). 


i, -iy +i, cosy —i, siny 
p =i, +i, cosy – p)—i, sin(y- e). (5) 
| =) ti, cos(y + p)-i, sin(y + p) 


-. 


~ 


1 
Wo aN ty, t y.) 


2 
Va 7a lv, cosy +y, cosly — p)- w. cosly + p)]; (6) 
2p | ; 
Va73 [va siny +y, siny — о) v. sin(y + p) 


where p= =120° for a machine with 


symmetrically arranged three phase windings, 


у=т+@ E — angle between the fixed phase axis a and 


the rotating longitudinal axis of the rotor. 


After the transformation, we obtain expressions for 
currents i,,i,,i, through flux linkages and trigonometric 


c 


functions of angle у: 
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і = (и, у, ++] 24-2 (a~ мат |+ 
31; 3. 3 


1 1 PENES | 
TW, -34* (а &)sin y * ^a - g)cosy siny = 


1 1 2. NS | 
DE Е 


— y ,b:cosy -WY aC cosy t v, h-siny 


(7) 
i bod О, 
i =z eve) Таз ila- одун y* 
3x, з 3 
V3 | 2 1 
+*=(a-g)cosy-siny TV, Za- -(a- g)4 
3 3 2 
1 v3 
*ala- s g)sin? ут 85 8 )соѕу :віпу 
ЕЕЕ n 
1 1 
ТЕ siny pr Е siny pr 
h 5 dn 438 as. 
Vahl іду + cosy 
(8) 


i. геу. ty, ни) еи, заз ila- eji y- 
Xo 3 3 


3 
V3 | 1 1 
- (а g)eosy-siny +0, ваа 


2 


ОСЕ > (a s)+=(a g)sin? y + 


v3 1 М3; 
+a- g)cosy siny ty b: ас + 


+ Cc: | cosy + 3 sin h | sin 1 COS 
V ua 2 у 2 Y TY kq 2 у 2 y 
(9) 


2 1 

i, 2——y,b:cosy +y,b 1cosy siny + 
3 3 3 

10) 


1 43. 
ty.b А * dy , = еу 


5 1 J3 
lico POS Pe сы йт: + 
(11) 


nee d ew, + fV, 
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i, E nh СЕ Ж + 
* og 3 3 
(12) 
1. X3 
ty hl QUEE ee + Куа 


The magnitude of the electromagnetic torque in the а, 
b, c axes is obtained after substitution into the well-known 
formula 


M, = V gla m» 
currents and flux linkages i,,/,,w,.V, expressed 


through i,.i, i, v ,, v ,, v, in the following form: 


243 5 >; jum 
м, =? Gai eu oL) -L0 


Substituting the expressions for currents i,,i, ,i, 


through flux linkages y,,w,,y,, we obtain the final 


a? 
expression for the electromagnetic torque used in 
calculating transient modes at the CT: 


243 1 
М, = By, «) (v, ty, +y,)+ 


3x, 


tV, $a-=(a~g)sin’ jp 


[1 1 ‚о, МЗ ; 
+] - a sa g)sin* y-+~=(a-g)eosy-siny |+ 


[11 so, МЗ | 
ty, –-а+ (а g)sin y — (а — в)соѕу:віпу – 
| 3 3 3 
-y jb-cosy ас cosy +,,h-siny}+ 


26 


1 
9 V. V, ) 3x, (Ж * V, + V. )+ 


1 1 TEE, | 
ЛЕЕ Sec enr 


2 1 1 А 
TE (a e) Ца g)sin? y- 
EC zi Lisci g) 
3 xi 2 


2 . 2 ( B Р 1 
- 3a - g)sin y |-W pb t+ Vc s siny ЕА 


1. vB 
Wight P 5 cosy |r + 
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2З | 


9 


1 
V. v) (v, +W, ty.) 
Xo 
ТЕЕ g)sin? yid - sony sin |+ 
3 3 3 
ty] dae g)-Za—g)sin® | 


2. 


DE y+ 


3 


ES h КА V3 oos 
2 Ай б И? 2 7 2 "arc 


Ш. TAKING INTO ACCOUNT THE EXTERNAL 
REACTIVITY OF THE SYSTEM 


NEJ | 1 
+—(а—в)созу-вїпу + ,Ь+уыс aad 


The system of equations (1) is compiled for the 
operation of a synchronous machine in parallel with an 
infinite power system. In this case, e,,e,,e,- are the 


a? "b? 


voltages on the tires of infinite power. 


The reactance and active resistances of transformers 
and lines from a synchronous machine to infinite power 
buses are included in the leakage reactance of the stator 
winding and in the active resistances of phase windings. If 
an accident occurs somewhere in the middle of a power 
line, then the voltage at the place of the accident will differ 
from the voltage of the buses of infinite power. In this 
case, part of the reactivity of the line from the generator to 
the point of the short circuit is added to the leakage 
reactivity and the machine is connected to the network 
Z,,Z3,Zc, which 


through additional resistance 


represents the resistance of the remaining part of the line 
from the point of the short circuit to the infinite power 
buses with voltages e,,e,,e- that do not depend on the 


mode. Let us consider a special case when z(p)=r+.xp. 
The system of voltages e,.e,,e, applied to the phases of 


the stator winding is determined from the Kirchhoff 
equations and depends on the mode under consideration. 
Kirchhoffs equations with the assumption in Fig. 1 emf 
direction and currents have the following form: 


€, = ед 1424 = ед 14ГА — XA la 
e, = —6g - lgZg = ер — lgl'g — Xgplg p. (14) 
€, — ес 14 = ес Z lere T Xe plc 
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Fig. 1: The direction of emf and currents in system 


After substituting (14) into (1), we obtain equations for 
the analysis of asymmetric modes, taking into account the 
resistances between the generator and the system. 


Phase voltages on buses of infinite power are 
determined by the well-known expressions: 


e, —e:Sinc 


cose . (15) 


COS | 


IV. ADDITIONAL RELATIONS FOR 
CALCULATING SOME ASYMMETRIC MODES 


А 1. 
ey = e:sin(z p)- (o 


о! у SA 


; 12. 
ec =e-sin(t + p)= E 


To solve system (1) taking into account (14), it is 
necessary to establish a connection between the currents of 
the generator i,,i,.i, and the currents in the buses 


i,,ip,ic or introduce additional relationships 


characterizing the operating mode of the synchronous 
machine. 

1. Asymmetry of the resistance system between the 
generator and the system Z, #Z, #Z,. In this case, as 
follows from Fig. 1 

i, =ty3, =1g5t, T dc. (16) 

System of equations (1) taking into account (14) is 

simplified 
PY, = ед - i.(r, + r,)- x4 pi, 
PW, =e, ~i,(x, + r,)- Xy pi; .(17) 
ру. =-—ес -i (r, +ro)-xepi, 
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Fig 2: Single-phase generator phase short circuit to 


common neutral generator-system 


2. Single-phase short circuit of the generator phase to 
the common neutral generator-system. In general, the short 
circuit mode can be calculated in the presence of resistance 
asymmetry Z, £ Zg £ 2с. 


Additional relations are obtained from Fig. 2. 


i, 2 i5 ;i, - ic ;e, =0. (18) 


The system of equations (1) takes the form: 


pV, = m^ 
pv, = -e, —i,(r, +) xs pi, . (19) 
pV,--ec - i (r, + ғ.) Xc pi, 
Short circuit of two phases of the generator to a 
common neutral generator-system. As follows from Fig. 3, 


when phases a and b of the generator are shorted to 
neutral, we obtain the following additional relations: 


i, ic ;e, = 0;е, =0. (20) 


Fig. 3: Short circuit of two phases of the generator to a 


common neutral generator-system 


Taking them into account, system of equations (1) 
takes the form: 


py, = cir, 
PY, = yr, . (21) 
PY. = ес ~i,(r, + rc )- xc pi, 

4. Two-phase short circuit. With a short circuit of 


phases a and b, we obtain only one obvious relationship 
from Fig. 4: 
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Fig. 4: Two phase short circuit 


і =. (22) 


с 


The other two relations are obtained after introducing 
an additional assumption: 


The values of active and reactive resistances between 
the short circuit point and infinite power buses in short- 
circuited phases are equal to 


Z,7Zgy = Рс, ТА 9 Tg o lc ХА — Xg = Хо. 


The introduction of this assumption allows us to obtain 
differential equations for stator circuits, which are easily 
solved by numerical methods. 


From expressions 


€, — €, б: | (23) 


e, — ёр — IBZ a 


we obtain expressions for e, +е,и e, —e,: 


€, te, --(e, вез), +1) 
: (24) 
е„—е„=0 
Using known relations e, +ев t ec = Оапа 
i, ig tic =0, get 
€, +e, -—ec Hi, tig = і. (25) 


After substituting (25) into (24) we get 


€, t€, = ес + Z4; 


e, е = 0. 
From here 

е Z 

e, = С +241; 
2 2 
е Z 

e, C qp A. Е 
2 2 


The system of equations (1) takes the form: 


ес X AI. , A 
= + 0r c s 
PY, 2 2 | di 2 
€ х, рі LY, Е 
py, = C AP „Са irs) (26) 


2 2 
DV, — ec = Xc pi, -i, r +). 
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5. Phase failure a. 


Phase loss is a special case of asymmetry of the 
resistance system, when the resistance 15 „=. 


However, it is not possible to introduce such a resistance 
value into the CT. We can limit ourselves to including a 
finite, but sufficiently large value Z, in the phase. Its 


value is limited by the scale factors adopted for 
resistances. 


This technique can also be used to consider the 
operating mode during phase failure. Consider the current 
in phase a to be equal to zero i, =i, =0, which 


corresponds to the physics of the process, and the flux 
linkages of phase a, determined by solving the first 
equation of system (1), are found using the equation 


V, Е IX ap ЕЗ 1. Хас + Dp Xy + аа * Ua X akg , 
where Xap» Xac>Xaf > Xazg — mutual reactivity of phase а 


and the corresponding circuits, determined by the 
following relations through the known parameters of the 
machine: 


х -in (а +X ix. - X, Jeos(2y - р) 


1 1 1 
„=д®-(х,+х,)-Щх, ox esr р). en 
XQ = Xaa * COSY 
Хак = Хад siny 


6. Single-phase short circuit of the generator to the 
neutral of the system (Fig. 5) is accepted: r, =r, 


a р 


= Fo , 
locked phase а. 
We obtain relations for currents 


i, =i, i, =i¢ (28) 


Fig. 5: Phase failure a 


To obtain the third necessary relationship, we compose 
the Kirchhoff equations for stresses: 


Tir ir, —€, —1,Z, = 0; 
PY, aa PY, › b ёв a B (29) 
PY, + lafa PY. Ut, ес 1.4с < 0. 


Adding equations (29) and using the relations: 


ру. +руь+ру. = 0; i, +i, +i, = 0. (30) 
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get: 
Зру, + 3igh, = ес +і Zc teg +12, 


where 
1 н А ; 
PW = rigs tide +е» 12) - ir (31) 


The equations for phases b and c are obtained 
similarly: 


1 Р . | 
PY, = -30e, -eç +2i,Zg -izc)- irs 
(32) 


1 | ; ; 
РҮ re -zlec — epg + 21.2 =). 


7. We consider a two-phase short circuit of the 
generator to the system neutral (Fig. 6) to be equal to the 
active resistance of the generator phases r,-7r,-r,; 


с 


Phases a and b are shorted. 


Fig. 6: Single-phase short circuit of the generator to neutral 


Composing, as in the previous case, the Kirchhoff 
equations 


Fig. 8: Two-phase short circuit to generator neutral 


DV, Hiza ру. ir. ec Ж =0 
pV,tir-py,-in-0 (33) 
DV, tir,- py, ir, -ec іс = 0 


get 


8. Single-phase short circuit to generator neutral 
(Fig. 7). Phase a is short-circuited. For phases b and c with 
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a ratio of i, =i, and i, = і. , from the Kirchhoff equations 


we obtain: 


PY, = mL 
DV, = ед +1424 ер —lyZg lh 
PY. ES ед +1424 ес lic 1.1, 


ід = -(i, ic) 


. (35) 


Fig 7: Single-phase short circuit to generator 


9. Two-phase short circuit to the neutral of the 
generator (Fig 8). Phases a and b are closed. For phase c, 
the relation holds i,=i,. Compose the Kirchhoff 


equations (Wood & Wollenberg, 1996) (Yakimchuk, 
2000) (Yan & Xu, 2020). 


ру, =i, 
DV, = iy, 
i, +i, ti, =0 . (36) 


РҮ tir tecticzo-e,—-i,z, =0 


PY. tir. tec ticzc — ez —igZg =0 


From the last two equations, under the condition 
Z, =Z, = Z, , we obtain 


3 | ; 
PY. = -5 lec іол), (37) 


Taking into account the reactivity between the 
generator and the infinite power system led to the 
appearance of derivatives of the stator phase currents on 
the right sides of the differential equations of phase 
voltages. 


V. SIMULATION RESULTS OF ASYMMETRIC 
TRANSIENTS 


The results of analytical studies of the considered 
asymmetrical transient processes (with asymmetrical short 
circuits) on phase coordinates are summarized in Table 1. 
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Table 2 and Fig. 9 present the numerical results of 
computer calculations based on the obtained solutions for 
transient asymmetrical short-circuit processes in the 
electrical system. The table data confirms the effectiveness 
of the proposed approach to studying asymmetric transient 
processes on phase coordinates. 


= ]U, р.а. 


13572911131517192123252729313335373941 434547 49 


Fig. 9: Voltage profile for a synchronous generator- 
grid network with asymmetrical short circuits 


As can be seen from Table 2 and Fig. 9, in case of short 
circuits on buses 28, 32 and 39, the voltage on these buses 


U =0 and in consumer nodes 1-9, 25-32 and 36-50 active 


and reactive loads is zero. On buses 40-47 the voltage is 
within (0,520 —0,594)U The results 


nom 


obtained 


confirmed the pre-modeled modes for the asymmetric 
transient modes considered above 


International Journal of Advanced Engineering Research and Science, 10(11)-2023 


VI. CONCLUSION 


1. Using phase coordinates, a method for calculating 
asymmetrical transient modes of a synchronous generator 
operating on an electrical network is presented. As a result 
of transformations, differential equations of a synchronous 
machine were obtained in axes a, b, c for stator quantities 
and axes d, q, O for rotor quantities. The variable 
coefficients included in these equations are functions of 
siny and cosy , therefore their calculation on a CT does 


not greatly complicate the calculation task compared to the 
calculation of symmetric modes in the d, q, О axes. 


2. Equations are given that make it possible to take into 
account, under certain assumptions, active and inductive 
resistance, the connection between the generator and the 
system of infinite power. Based on the conditions 
characterizing asymmetrical modes, equations for the 
stator circuits of a generator operating on infinite power 
systems for various asymmetrical short circuits are 
obtained. 


3. Numerical results of computer calculations based on 
the obtained solutions for transient asymmetrical short 
circuit processes in the electrical system are presented. The 
table data confirms the effectiveness of the proposed 
approach to studying asymmetric transient processes on 
phase coordinates. 


Table 1: Results of analytical studies of asymmetrical transient processes 


Types of asymmetrical 


System equation (1) 


short circuits Я 


PY, = e, — 124 =, PY, = €, – 1,25 n PY, = ec -i.e E 
Single-phase short circuit to py, =L, pV,-—-—eg-ir,—iyzg | DV, ——ec ісі, 
common neutral 
Two-phase short circuit to pv,--ir, PY, = і, PY. = Cc -i,zo - i,r, 
common neutral 
Two-phase short circuit e iz ec dz ру. = есі і 
С . A C A . c C c*A 
PY, tne PY, = +t- 17 = 
2 2 2 
Single-phase short circuit to 1 1 1 
ру = (6. + с + py, = (е6 -ec + ру = (06-е, + 
system neutral ЭЕ: d 3 d 3 


+ 2i,2%¢ – 1,25 )- ir, 


+ 2i, Zp ш їс )- int, 


Two-phase short circuit to 
system neutral 


1 2 
PY, = zlec +i zc) ib, ру, =-Flec ti.zc)- ir, 


Single-phase short circuit to 
generator neutral 


DV, —e,*i,z,—ec— 


ісік, = ед + 


с 


PW, =ед +1424 
ёр 14р yl, 


+1524 ес - i (z, * zc)- ir, 
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Two-phase short circuit to 
generator neutral 


PY, ==, 


PY, =n, 


3 ; : 
ру. = -F lec ticzc)- ir. 


Table 2: Calculation results based on the obtained solutions for asymmetrical short circuit 


Node | Node U Angle P, Qe Pi Qi 
No. type | (p.u.) (deg) (MW) | (MVAR) | (MW) | (MBAP) 
1 3 1,050 -1,916 0,000 0,000 0,000 0,000 
2 3 1,051 | -121,674 0,000 0,000 0,000 0,000 
3 3 1,055 116,154 0,000 0,000 0,000 0,000 
4 3 1,052 68,661 0,000 0,000 0,000 0,000 
5 3 1,050 -30,575 0,000 0,000 0,000 0,000 
6 3 1,062 | -150,839 0,000 0,000 0,000 0,000 
7 3 1,053 87,380 0,000 0,000 0,000 0,000 
8 3 1,056 -32,461 0,000 0,000 0,000 0,000 
9 3 1,056 | -152,688 0,000 0,000 0,000 0,000 
10 6 1,046 87,388 0,000 0,000 9,111 5,466 
11 6 1,044 -33,011 0,000 0,000 9,089 5,453 
12 6 1,037 | -152,765 0,000 0,000 8,968 5,361 
13 6 1,041 87,096 0,000 0,000 10,833 6,861 
14 6 1,039 -33,436 0,000 0,000 10,791 6,834 
15 6 1,029 | -153,101 0,000 0,000 10,587 6,705 
16 6 1,042 87,080 0,000 0,000 12,676 7,965 
17 6 1,045 -32,377 0,000 0,000 12,625 7,934 
18 6 1,032 | -153,073 0,000 0,000 12,428 7,809 
19 6 1,040 87,025 0,000 0,000 10,820 6,352 
20 6 1,038 -33,469 0,000 0,000 10,778 6,525 
21 6 1,029 | -153,152 0,000 0,000 10,588 6,706 
22 6 1,040 87,014 0,000 0,000 14,420 9,012 
23 6 1,038 -33,458 0,000 0,000 14,368 8,980 
24 6 1,029 | -153,148 0,000 0,000 14,129 8,830 
25 1 1,070 0,000 38,839 23,199 0,000 0,000 
26 1 1,070 | -120,072 32,511 21,553 0,000 0,000 
Table 2 (continue) 
Node | Node U Angle Р, о, Pi Qi 
No. type | (ри) (deg) (MW) | (MVAR) | (MW) | (MBAP) 
27 1 1,070 | 120,072 | 37,072 16,920 0,000 0,000 
28 5 0,000 0,000 0,000 0,000 0,000 0,000 
29 4 1,070 | 90,856 | 24,543 11,722 0,000 0,000 
30 4 1,070 | -29,168 | 16,229 13,057 0,000 0,000 
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31 4 1,070 | -149,192 | 19,228 5,183 0,000 0,000 
32 5 0,000 0,000 0,000 0,000 0,000 0,000 
33 6 1,055 | 87,412 0,000 0,000 6,680 2,598 
34 6 1,058 | -32,443 0,000 0,000 6,710 2,610 
35 6 1,058 | -152,634 | 0,000 0,000 6,714 2,611 
36 4 1,070 | 85,282 10,083 9,429 0,000 0,000 
37 4 1,070 | -31,670 | 8,970 7,914 0,000 0,000 
38 4 1,070 | -151,694 | 10,947 7,732 0,000 0,000 
39 5 0,000 0,000 0,000 0,000 0,000 0,000 
40 3 0,520 | -93,147 0,000 0,000 0,000 0,000 
41 3 0,863 | 82,090 0,000 0,000 0,000 0,000 
42 3 0,751 | -9,572 0,000 0,000 0,000 0,000 
43 3 0,192 | -72,025 0,000 0,000 0,000 0,000 
44 3 0,297 | -167,711 | 0,000 0,000 0,000 0,000 
45 3 0,601 | 67,024 0,000 0,000 0,000 0,000 
46 3 0,599 | -33,474 | 0,000 0,000 0,000 0,000 
47 3 0,594 | -153,152 | 0,000 0,000 0,000 0,000 
48 3 1,034 | -122,883 | 0,000 0,000 0,000 0,000 
49 3 1,042 | 116,762 | 0,000 0,000 0,000 0,000 
50 3 1,032 | -3,459 0,000 0,000 0,000 0,000 
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